Quantum criticality in a Mott pn-junction in an armchair carbon nanotube 
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In an armchair carbon nanotube pn junction the p- and n- regions are separated by a region 
of a Mott insulator, which can backscatter electrons only in pairs. We predict a quantum-critical 
behavior in such a pn junction. Depending on the junction's built-in electric field E, its conductance 
G scales either to zero or to the ideal value G — 4e 2 //i as the temperature T is lowered. The two 
types of the G(T) dependence indicate the existence, at some special value of E, of an intermediate 
quantum critical point with a finite conductance G < 4e 2 /h. This makes the pn junction drastically 
different from a simple barrier in a Luttinger liquid. 

PACS numbers: 71.10.Pm, 64.70.Tg 



Transport measurements in carbon nanotube (CNT) 
devices reveal strongly correlated behavior of conduc- 
tion electrons. Experiments on tunneling into single-wall 
CNT [U 12 or across a barrier interrupting a CNT [2] 
demonstrated a power law bias- and temperature- de- 
pendence of the current, consistent with a gapless ex- 
citation spectrum and Luttinger correlations of the elec- 
tron liquid. Recent experiments [3] showed that armchair 
CNTs develop a gap in the spectrum of charge excitations 
at zero doping. The gap formation could be attributed 
to the electron Umklapp processes [H [5], which back- 
scatter pairs of electrons. These processes drive armchair 
CNTs into the Mott-insulating state. In contrast, un- 
doped semiconducting CNTs are simple band insulators. 
Although Mott and band insulators are qualitatively dif- 
ferent, no clear experimental signature of that dichotomy 
has been observed in CNTs. 

In this Letter we predict that the difference should 
manifest itself in the conduction of a pn junction formed 
by bipolar doping of a CNT. In contrast to a simple 
"band" pn junction, in which the barrier between the p- 
and n- regions is formed by a band insulator, a "Mott" 
pn junction formed in a pristine armchair CNT does 
not back-scatter single electrons. However, the Umk- 
lapp backscattering of pairs of electrons remains effective 
near its center. In short junctions these processes are ir- 
relevant and do not alter the perfect zero-temperature 
conductance (G = Ae 2 /h). In longer junctions, Umk- 
lapp processes bring about the Mott-insulator state and 
drive the zero-temperature conductance to zero. The 
bipolar doping of a CNT can be achieved by gating the 
nanotube [5J, which enables one to change the junction 
length by controlling the built-in electric field E. By 
varying E at a fixed low temperature T, one can tune the 
conductance G of the junction between almost perfect 
and zero values. Lowering the temperature makes this 
crossover sharper and culminates in a zero-temperature 
quantum phase transition occurring at some critical field 
value E — E* . The described dependence on E and T 
the Mott pn junction conductance, G(T, E), is drastically 
different from that of a band pn junction, where G van- 




FIG. 1: (a) Electron spectrum near the Dirac points a = ±1 
(+ or — indicates parity), (b) Near x = the gate potential is 
U(x) ~ —eEx and saturates to ±Uq in the p- and n- regions, 
(c) Schematic picture of the conductance dependence on E at 
different temperatures. The G(E) curves at different temper- 
atures intersect at the fixed point and increase monotonically 
with E. The step in G(E) at T = 0, indicates a quantum 
phase transition, (d) RG flow for the conductance. 



ishes at any E in the low temperature limit, T — > 0. We 
make detailed predictions for the G(T, E) dependence in 
the vicinity of the quantum phase transition. 

We consider a single pn-junction that is formed by 
a static potential U(x) imposed by two gates, as in 
Fig.jljb). In the p- and n- regions U(x) saturates to ±J7q, 
while near the center of the pn-junction U(x) ~ eEx, 
where — e is the electron charge and E is the built-in 
electric field. If the gate potential is not too strong, 
\Uq\ <C vf/R, where R is the CNT radius (we set h = 1 in 
the intermediate formulas), the pn-junction may be mod- 
eled using the low energy CNT band structure shown 
in Fig. [I] (a). Two spin degenerate bands intersect at 
two Dirac points (valleys, a = ±1). The electron wave 
functions in the intersecting bands have opposite parity 
P = ±1. They are either symmetric (P = 1) or antisym- 
metric (P = —1) in the AB sublattices [7j. We assume 
that U(x) is smooth on the interatomic scale and does 
not break the symmetry between the A and B sublat- 
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tices. Then both inter- and intra-valley backscattering is 
absent, and we may write the noninteracting part of the 
electron Hamiltonian as 



H a = 



E 



dxip+ r(T (x) [-irv F d x + U(x)]ip ara (x). 



Here vp ~ 8 * 10 5 m/s is the Fermi velocity, a — ± is the 
valley index, r = ±1 labels left and right movers, and a 
the electron spin. 

We assume that the number of unit cells N around the 
circumference of an (N, N) armchair CNT is large. To 
zeroth order in 1/N the Coulomb interaction of electron- 
electron (e-e) backscattering is absent [81. The forward 
scattering part of the e-e interactions is, 



V f = 



V(0) 



dxn 2 (x); n(x) 



^tra{ x )^<xra{x). 



Here V(0) ~ 2e 2 \n(d/R) is the forward scattering ma- 
trix element (we assume that the Coulomb interaction is 
screened by the gate at a distance d from the CNT). 

The gate potential U(x) results in a position- 
dependent doping density characterized by the Fermi 
wave vector kp(x) ~ —U(x)/vp. We can bosonize the 
electron operators by the standard procedure [9] 



Ipara (^) 



F, 



e ir U" k F (x')dx' 



* «(x)]+i6 Q „(i) 



(1) 



Here £ ~ R is the short distance cutoff, F aa are 
Klein factors. The bosonic fields <E> and 9 obey the 
standard commutation relations [$ a n{x),®a'cr'{x')\ — 
-iir5 aa >5 aa >9(x — x'). 



In the bosonic representation the electron density is 
n(x) — Akp(x)/ir — ^2 atT d x <^ arj {x)/iT and the forward 
scattering part of the Hamiltonian Hj = Hq + V/ reduces 
to the same form as in a uniform CNTs [TU] 



H f = J ^E^[^(W) 2 + (W 2 /^] 



(2) 



Here j — c±, s± labels charge (c) and spin (s) modes 
that are symmetric (+) or antisymmetric (— ) in the 
valley index a. They are related to the fields $ Qcr by 
®aa = [ c &c+ + a$ c - + v®s+ + acr$ s -] /2. Only the c+ 
mode carries charge and the other three are neutral. The 
mode velocities are Uj = Vf/Kj, where the Luttinger pa- 
rameters are K c+ = 1/ yl + 4:V(0)/nvp -C 1 and Kj = 1 
for the neutral modes. 

Backscattering interactions are small in 1/N [S], and 
may be treated as perturbations to the Hamiltonian 
Some of them are relevant and qualitatively change the 
low energy physics. The most relevant backscattering 
interaction corresponds to the so-called Umklapp pro- 
cesses, which scatter two right-movers into left-moving 



states or vice- versa [H [TU] ■ In the presence of doping the 
Umklapp Hamiltonian can be written as [5] 



H, 



dx 



cos 2<S> c+ {x) - 4 / k F {y)dy x 



{53 cos[29 s _ (x)} + (g 3 - gi ) cos[2$ s+ (x)] 

+ gi(cos[2$ c _(a;)] - cos[2$ s _(a;)])} . (3) 

Here the coupling constants 171 and 33 are of order e 2 /N 
and 53 > gi . The low energy electron Hamiltonian of the 
pn-j unction is given by the sum of Eqs. ^ and ([3|, 



H — Hf + Hjj. 



(4) 



The position-dependent Fermi vector kp(x) in Eq. ([3]) 
saturates to constant values ±feo in the n- and p- regions 
whereas near x = it has a linear dependence of x, 



kp(x) = —x/L 2 E . 



(5) 



Here the length scale Lp is defined by the built-in electric 
field E controlled by the gate voltages, Lp ~ y/vp/eE. 

Relating the current operator to charge field as j = 
—2edt§c+(x = 0, t)/"K and using the Kubo formula we 
can express the device conductance as 



8e 2 

G = i— limwg w (0,0). 

7m co-to 



(6) 



Here Qu(x, x') is the retarded Green's 
function of the charge field, Qu(x,x') = 
-i f™ die™* {[$ c+ (x,t),$ c+ (x',0)}}. 

In the absence of Umklapp interactions evaluation of 
the Green's function is straightforward and gives 

G°(x-x') = -i^(K c+ /2uo + )e l ^ x - x '\' u ^, (7) 

where u> + = oj + ir\. When substituting this expression 
into Eq. (|6| it should be born in mind that the dc con- 
ductance is controlled by the leads, where LL effects are 
absent. Thus the Luttinger parameter should be set to 
unity, K c+ -> 1 [UJ, yielding G — 4e 2 //i. 

The Umklapp processes degrade electric current by 
backscattering pairs of electrons and thus decrease the 
device conductance. They are most effective near the 
zero doping point x = and strongly suppressed deep in 
the p- and n- regions. At low temperatures, T vpka, 
the Umklapp backscattering in the p- and n- regions re- 
sult in exponentially small ~ exp(— 2vpka/T) resistivity. 
We assume that the length of the p- and n- regions is 
not sufficient to compensate for this exponential small- 
ness and neglect this contribution. In this approximation 
backscattering arises from the spatial region |ac| < T/eE, 
where |fc_p(x)| < T/vp. At T <C Uo the spatial depen- 
dence of kp{x) in this region is linear. Therefore the pn- 
junction may be modeled by the Hamiltonian Q with 
kp(x) given by Eq. |5| in the entire space. 
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Spatial localization of backscattering enables us to ex- 
press Green's function of the charge mode as 



can only scatter pairs of electrons. The Umklapp scatter- 
ing may be described by an effective "impurity" Hamil- 
tonian H 2 ~ Lp%u{Q). Upon reduction of the energy 



/oo " " 1 ' 1 

dxdx'T u (x,x')G°(0). (8) band width > A o ~ "Fit -> ^ A = A e-', new operators 
-oo arc exonerated in the effective Hamiltonian of the bar- 



Here T u (x, x') is the part of the T-matrix that corre- 
sponds to scattering of a plasmon into a single plasmon. 
Using Wick's theorem (c.f. Appendix A of Ref. [T^]) one 
can show that it can be expressed as 

/>oo 

T u (x,x') = -i dte iuJt ([d^n u (x,t),d<s,n u (x , ,0)}) 



(9) 



+ 5{x-x'){din u { x )) 



where Tlu(x,t) is the Umklapp Hamiltonian density and 
9$ denotes a partial derivative with respect to <fr c +- 

In the low frequency limit needed for evaluating the 
dc conductance the unperturbed Green's functions G® in 
Eq. (TSJ) are given by Eq. Q with K c+ = 1. Then the 
deviation 6G of the device conductance from the ideal 
value, 4e 2 //i, may be expressed as 



SG 



e 2 2n 
-i— hm — 
h w->o uj 



dxdx'T u {x,x'). (10) 



The T-matrix properties are dominated by the fluctu- 
ations of bosonic modes with frequencies on the order 
of the temperature T and characteristic spatial scale of 
Lt ~ vp/T. Provided the device length is longer than 
Lt the T-matrix needs to be evaluated using the Hamil- 
tonian Q with the values of the Luttinger parameters 
corresponding to the device interior. 

The energy gap A ~ (vf/£,)(93/ v f) 1 ~ K ' !+ [3] induced 
by the Umklapp interaction in a uniform CNT at zero 
doping defines an additional characteristic length scale, 
C = Vp/A. Backscattering at the pn-junction is weak for 
and strong for Lp >• £. 
For Lp <C C the correction to the ideal conductance 
may be expanded in perturbation series in Hu, Eq. fcfy. 
Using Eqs. (10 1, dOJ) we get to second order in Hjj , SG^ 2 ' = 
' c+ {L E /^ 2 {3g 2 -2g m + 2g 2 )/v 2 F , where 
i c+ + 1,K C+ + l)/2n, with B(x,y) being 



-c{2irTd/u c+ Y' 
c = (e 2 /h)B(K c 
the Euler Beta function. This result may be rewritten as 



6G^ 



2 /h)(L E /0 2 (T/A) 



1K a 



It vanishes at T 0, 



which corresponds to irrelevance of point-like Umklapp 
scattering in the renormalization group (RG) sense. 

Although higher order terms in perturbation theory 
are smaller in powers of LpjC, they have different tem- 
perature dependence. In the fourth order correction to 
the ideal conductance the most relevant term is SG^ ~ 
-{e 2 /h)(L E /0\T/Af K ^- 2 . 

This result can also be obtained using the following RG 
considerations. Because of localization of backscattering 
to the the region of size |x| < Lp at low temperatures, 
T <C vp/Lp, the pn-junction acts as a point scatterer. In 
contrast to a potential point scatterer [13] the pn-junction 



are generated in the effective Hamiltonian of the bar- 
rier. They have the form of higher powers of H^. The 
most relevant of those is H 4 ~ Ar 4 cos (4$ c+ (0)), with 
the scaling dimension AK C+ — 1. The correction SG^ 
can be obtained by the lowest order perturbation theory 
in H 4 , 5G& ~ (e 2 /h)jtdt([H 4 (t),H 4 (0)}) cx T 8 ^+~ 2 . 

For K c+ < 1/4 the H 4 is relevant and SG^ diverges 
at zero temperature. The RG fixed point (FP) at perfect 
transmission is unstable, and the system flows to strong 
backscattering even at Lp <C C We show below that the 
strong backscattering FP is stable. The expected RG 
flow is as shown in Fig.JTJd). It corresponds to vanishing 
zero temperature conductance. We note that the tem- 
perature dependence of conductance is nonmonotonic. In 
the regime of applicability of perturbation theory 5G = 
-a (Le/C) 2 (T/A) 2i ^+ - b (Lp/Cf (T/A) 8i ^+- 2 , where 
a and b are constants on the order of e 2 jh. The maximum 
conductance is reached at T ~ A (L p / Q) 1 ^ 1 ~ 5Kc+ \ 

Because of the large value of the interaction constant, 
e 2 jvp sa 2.7 and the sensitivity of the forward scattering 
matrix element V(0) to screening by the gates both cases 
K c+ < 1/4 and K c+ > 1/4 may be realized. 

For K c+ > 1/4 all backscattering operators generated 
in the process of renormalization are irrelevant. There- 
fore the FP at perfect transmission is stable for weak 
backscattering, < (. At the same time, the anal- 
ysis below shows that a perfect reflection fixed point 
is stable at K c+ < 1. This indicates the existence at 
1/4 < K c+ < 1 of a quantum phase transition in the con- 
ductance of pn-junction, controlled by the ratio Lp/C,. 

At strong backscattering, Lp >• (, the pn-junction 
may be viewed as two semi-infinite Luttinger liquids sep- 
arated by a strong barrier where the bosonic fields are 
pinned to one of the classical minima of the Umklapp 
Hamiltonian. The minima form a periodic lattice in the 
four-dimensional space $>j = 7rn 3 /2, where rij are inte- 
gers, which are either all even or all odd. Charge trans- 
port between the p- and n- regions proceeds via tunneling 
between different points in the lattice [13] • The tunneling 
operator corresponding to the shift {nj} — > {rii + Srij} is 
exp{i Y2j[0j(xi) — 9j(x2)]Srij /2]}, where x\ and x-i label 
the points just to the left and to the right of the barrier. 
Its scaling dimension is 1 — J2ji^ n j) 2 /i^j)- The oper- 
ators that transfer charge across the junction must have 
nonzero 6n c+ . For K c+ < 1 all of them are irrelevant, 
and thus the zero temperature conductance vanishes. 

The stable FPs at perfect transmission and reflection 
for 1/4 < K c+ < 1 must be separated by an unstable 
FP with an intermediate value of zero temperature con- 
ductance. The pn-junction can be tuned to this FP by 
adjusting the parameter Lp. The dependence of con- 
ductance on Lp at low temperature is schematically pre- 
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sented in Fig. [ljc) . The critical conductance and critical 
exponents at the intermediate fixed point depend on the 
values of the Luttinger parameters Kj . 

To analyze the intermediate FP near K c + — 1/4 we 
note that only two of the operators generated in the 
course of renormalization have scaling dimensions that 
are reasonably small: and Hjj (respective scaling di- 
mensions AK C+ — 1 and K c+ ). All other operators are 
strongly irrelevant. The proximity of the values of K c+ 
corresponding to marginality of Hjj and H 4 {K c+ = 0, 
and K c+ = 1 /4) enables one to use a kind of e-expansion 
in which both e = K c+ and £4 = K c+ — 1/4 are assumed 
small. In this case the intermediate FP is perturbatively 
accessible. Writing the effective pinning Hamiltonian 
as H eff = A[r 4 cos(4$ c+ ) + cos(2$ c+ ) Y^ a r a cos(2$ )], 
where r a are dimensionless amplitudes and a labels 
the neutral modes, we obtain the RG equations to lead- 
ing order in e, 64, 



dir a = -er a + r a r 4 , 



(11a) 
(lib) 



where d; = d/dl. The intermediate FP exists for £4 > 
and is located at r| = e and J2 a r *a = 16f4e- The FP con- 
ductance can be found by using lowest order perturbation 
theory with FP values of the reflection amplitudes [15] . 
G* = A(e 2 /h)(l - e 4 - e 2 /32e 4 ). The presence of the e 4 



in the denominator in the last term signals that Eq. ( 11 ) 
is valid provided 64 3> e 2 - Near the FP the conductance 
behaves as G(T) - G* ~ (e 2 /h){l - L E / 'L* E ){hv F /TQ X + 
where A + = — 2e4+2y / e 2 + 2ee 4 . A similar analysis in the 
limiting case of 1 — Kj — rj — ► 0+ gives the intermediate 
FP conductance G* oc rj 2 . Interpolation between these 
two limiting cases gives the RG flow shown in Fig. [ljd). 

The above picture is modified if single-electron 
backscattering is present. In a symmetric armchair CNT 
it can be caused by a magnetic field B applied along the 
CNT axis or by the electron-phonon (ep) interactions. 
The corresponding backscattering Hamiltonian may be 
written as [PoHTT] 

SH = -i dx y"yi riJ (x)tp a - ra r[aA u {x) + A B ], (12) 



where A B and A u are respectively the gaps induced 
in the single particle spectrum by the magnetic field, 
As = (n/2)eBv F R, and by the lattice deformation. At 
low temperatures only the twist acoustic (TA) phonons 
are important and A u (x) — gTd x u(x), where u(x) is the 
TA mode displacement and g F ~ 1/vN 1 is the cor- 
responding coupling constant. As in the case of two- 
particle processes, single electron backscattering is effec- 
tive only in the vicinity of the pn-junction, |x| < L B . 
The zero-transmission FP is stable with respect to single 
electron backscattering. 



For Le « ( and at small B and gr the stability 
of the FP at perfect transmission with respect to sin- 
gle particle backscattering may be examined by per- 
turbation theory. Bosonizing the fermion operators in 
Eq. ( [12] ) as in Eq. ([I]) we find the lowest order correc- 
tion SG ~ -(e 2 /h)[(L E /Z) 2 (g2/ps 3 T )(ZT/v F )( 3 + K ^ + 
{L E A B /v F ) 2 (v F T/^ K "+- i y 2 }. Here p is the mass per 
unit length of the CNT and st is the speed of sound for 
the TA mode [IS]. The temperature dependence of the 
two terms here shows the stability of the perfect trans- 
mission FP with respect to the ep backscattering and its 
instability with respect to B 7^ 0. 

Combining the stability analysis of the perfect trans- 
mission and perfect reflection FPs we conclude that the 
intermediate FP is not destroyed by the ep interactions. 
However application of a magnetic field drives the zero 
temperature conductance of the junction to zero. This 
reveals that magnetoresistance of a CNT pn-junction [18] 
is strongly enhanced by the Luttinger liquid effects. 

In conclusion, the Mott-insulating state of an armchair 
CNT should manifest itself in the dependence of the pn 
junction conductance G on the temperature T and built- 
in electric field E. In a broad range of parameters, the 
G(T,E) dependence is controlled by a T — quantum 
phase transition occurring at some value of E = E*. The 
critical behavior of G(T, E) is controlled by an unstable 
fixed point with conductance G* < 4e 2 /h. That suggests 
the possibility of the scaling analysis of experimental data 
used in investigations of quantum phase transitions [19] . 
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